The Distribution of Nearby Stars in Velocity Space 
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Abstract. From the tangential velocities of stars all over the sky, one 
can, in a statistical way, infer their 3D velocity distribution. An applica- 
tion to Hipparcos data reveals rich structure in the planar stellar motions: 
there are several moving groups, increasing in number with stellar type 
' but decreasing in importance. A distinct group of outward moving stars 

, with low rotation velocities might be associated with the Galactic bar: 

fSj ' growing a bar in a smooth model for the Galactic disk results in such a 

group, provided the Sun is outside the outer Lindblad resonance (OLR) 
and the orientation angle with the bar is in the first quadrant. 

The vertical motions show less structure, but the mean v z increases 
\ for large rotational velocities, which can be nicely explained by a warp of 

CO ■ the outer stellar disk. 
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1. Inferring f(v) from Hipparcos Data 
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^"'i ESA's Hipparcos mission measured the positions, parallaxes, and proper motions 

of many nearby stars. The radial velocity (RV) of most of these stars is unknown 
-£3 . (using only stars with known RVs gives a sample strongly biased towards high-?; 

stars [Binney et al., 1997]), and we cannot evaluate their space motions v, but 
only the tangential velocities Vf. However, using stars from all over the sky, we 
can infer their distribution f(v) in a statistical way, for instance, by maximizing 
^ ■ the log-likelihood 

£(/) = N- 1 Y^hxf dv r f( V = V t>i + V r h) (1) 
i 

with r i the unit vector in direction to star i. Since £(/) is unbounded, it 
must be regularized: Q a = C — aS is maximized instead with S(f) a penalty 
functional measuring the roughness of f(v) and a the smoothing parameter. The 
maximum of Q a (f) subject to / > and fdfivf=l is the maximum penalized 
likelihood estimate (MPLE). The degree of smoothness, determined by a, may 
be optimized by cross-validation using Monte-Carlo simulations (for more details 
and the numerics, see Dehnen, 1998, hereafter D98). 

This technique was applied to the kinematically unbiased sample of Hip- 
parcos stars defined by Dehnen Sz Binney (1998, see also Binney, this volume). 
I analysed seperately four colour bins of main-sequence (MS) stars, the non-MS 
stars, and the full sample. For the latter, optimal smoothing was used, while 
for the subsets, f(v) was slightly undersmoothed to preserve features that are 
likely to be real because they are present in more than one distinct subset. 
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Figure 1. v z measured for the full sample as function of v y . The A 
indicates the LSR. Lines represent models for the Galactic warp. 

Below I use a right-handed cartesian coordinate system with e x , e y , e z 
pointing towards (£, b) = (0°, 0°), (90°, 0°), and b = 90°. 

2. Vertical Motions: Evidence for a Warp in the Stellar Disk 

The projections of f(v) onto the v x v z or v y v z plane (Figs. 4&5 in D98) do 
not show any significant structure (see Binney, this volume, for reasons why). 
However, there is a skewness in the sense that v z increases for large positive v y , 
see Fig. |l[ This effect results naturally from a warp in the outer stellar disk: 
stars with large v y visit us from large R, where the warp causes them to move 
upwards. The lines in Fig. [l] are due to a simple quantitative model (see D98 
for details) for this effect. A warp starting within Rq (dashed) cannot fit the 
data, while warps starting at ~ Ro (solid) can. There is a degeneracy between 
the warp's amplitude and pattern speed f2 w in the sense that 

4 - 6(zi /kpc) 

kms-ikpc- 1 ~ 0.3(zi /kpc) 1 ' 

where z\q is the amplitude at R = lOkpc. For Hi, z\q ~0.3-0.4kpc, which would 
give Q w ~ 13 to 25kms -1 kpc~ 1 (retrograde w.r.t. rotation). 

3. Planar Motions: Moving Groups 

The projection of / onto the v x v y plane shows significant structure, see Fig. |2[ 
Surprisingly, there are, apart from the well-known moving groups, also several 
minor features, mainly at high v, i.e. on eccentric orbits, and for late stellar types 
only, i.e. made of old stars. The standard picture of moving groups forming via 
dissolution and mixing of stellar clusters cannot account for these features. A 
viable explanation, though, is a dynamical formation of the moving groups, by 
means of resonant locking. If the Galactic potential changes slowly, the resonant 
islands slowly sweep through phase space. Stars may then become locked to it 
and shifted to other phase-space regions (cf. Sridhar & Touma, 1996). In this 
picture, moving groups could be made of stars of different origin and age. 
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Figure 2. Distribution in planar velocities (v x : towards £ = 0°; 
v y : towards £ = 90°) for the four colour bins of MS stars (a-d), 
the non-MS stars (e), and the full sample (f). Velocities are rel- 
ative to the Sun, and the A indicates the LSR. Contours contain 
2, 6, 12, 21, 33, 50, 68, 80, 90, 95, 99, and 99.9% of all stars. Only 
the full sample (f) is optimally smoothed, the subsets are under- 
smoothed and contain artificial structures. However, features present 
in two or more subsets are likely to be real, e.g. the double peak at 
v x = — 25, — 45kms~ 1 and v y = — 50kms _1 present in (d) and (e). 
Note the four major peaks, best visible in the upper right panel, which 
correspond to (from left to right) the Hyades, Pleiades, Coma, and 
Sirius moving groups. At later stellar types there are also many minor 
moving groups, predominantly at high velocities (w.r.t. the LSR) i.e. 
on eccentric orbits, but the total number of stars associated with them 
diminishes with stellar type. 
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Figure 3. Top: as Fig. ^ but (d) 
and (e) analyzed together. Bot- 
tom: f(v) obtained from mod- 
eling: a bar has been grown in 
an initially smooth, axisymmet- 
ric, and exponential disk. The co- 
rotation radius is 4kpc, while solar 
radius and bar position angle are 
(i2o, = (7.5 kpc, 10°) for (b) and 
(8.0kpc,40°) for (c). 
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4. Planar Motions: Influence from the Galactic Bar? 

There is a distinct group of stars at v y ~ — 50 and —100 5a v x ^ 20 (in kms -1 ), 
see Fig. |3|a, yielding v x <0 at v y <— 30kms _1 , also found by Raboud et al. (1998) 
from a sample of Hipparcos stars with known RVs. These authors report high 
Z for these stars and argue they are a signature of the Galactic bar, as in Fux's 
(1997) models v x < at bar position angles <f> G (0°, 90°). Such a feature is created 
naturally by a bar grown in an initially axisymmetric disk, Fig. ||b&c, if we are 
outside the OLR. In this scenario, the outward moving stars visit us from inside 
the OLR, where the orbits move outwards for (f> S (0°, 90°). The position of this 
group in w-space depends sensitively on <f) and our distance to the OLR. 
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